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Phenomena  associated  with  acoustic  cavitation  in  liquids  are 
observed  to  change  abruptly  in  magnitude  with  small  changes  in 
the  amplitude  of  the  acoustic  pressure  or  intensity.  The  amplitudes 
at  which  changes  occur  are  called  cavitation  thresholds.  Investiga- 
tors have  noted  that,  once  such  thresholds  have  been  exceeded,  the 
phenomenon  being  observed  rises  to  a maximum  and  then  decreases. 

L.  D.  Rozenberg1  was  perhaps  the  first  to  point  out  the  existence  of 
these  maxima. 

Many  experimenters  have  determined  both  the  existence  and 
magnitude  of  cavitation  thresholds,  using  erosion  of  solids, 

cavitation  noise,  chemical  reactions  and  biological  effects  as  the 

. 2 
cavitation  phenomena  under  observation.  The  work  of  Kaufman, 

. 3 

Miller,  Griffiths,  Ciaravino  and  Carstensen  and  of  Clarke  and 

. 4 

Hill  may  be  cited  as  examples  of  the  use  of  biological  effects  m 

5 

threshold  measurements,  while  that  of  Neppiras  and  Coakley  may  be 
cited  as  an  example  of  use  of  cavitation  noise.  For  brevity,  the 
term  "cavitation  activity"  will  be  used  to  denote  the  wide  range  of 
phenomena  associated  with  cavitation  fields  or  zones  and  used  by 
various  experimenters  to  demonstrate  the  existence  of  thresholds  as 
functions  of  pressure  or  intensity. 

Attention  will  be  limited  here  to  the  remarkable  papers  of 

g 

Monakhov,  Peshkovskii,  Popovich,  Fomichev,  Chinyakov  and  Yakovlev 

. . 7 

and  Brandt,  Yakovlev  and  Peshkovskii  . These  workers  observed  two 
cavitation  thresholds  — one  for  the  onset  of  cavitation  activity 
and  a second  for  a marked  decrease  in  cavitation  activity  (except 
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for  the  radiation  of  noise) . Monakhov  et  al  distinguished  three 

regimes  of  cavitation.  Below  the  first  threshold,  there  was  no 

cavitation  activity.  Above  the  second  threshold,  cavitation 

activity  (as  evidenced  by  erosion,  for  example)  virtually  ceased, 

large  bubbles  were  observed  and  the  liquid  resembled  boiling  water. 

Monakhov  et  al  carried  out  their  experiments  at  a single  acoustic 

frequency  (17.8  kHz)  in  water  that  presumably  contained  gas  "seeds" 

or  "nuclei"  having  a wide  distribution  of  initial  radii.  In  the 

2 9 

terminology  suggested  by  this  author  ' the  first  regime  of  Monakhov 
et  al  would  be  a zone  containing  only  stable  cavities  while  the 
regime  bounded  by  the  two  thresholds  would  be  a zone  dominated  by 
transient  cavities.  In  the  third  regime,  the  zone  would  again  be 
dominated  by  bubbles  that  behave  much  as  stable  cavities. 

An  objective  of  this  paper  is  to  seek  a dynamical  basis  for 
the  existence  of  such  thresholds  and  of  the  observe^,  maxima  in 
cavitation  activity.  This  paper  is  the  third  in  a series  on 
cavitation  dynamics.  In  the  first8,  hereinafter  referred  to  as  CD: I 
a mathematical  formulation  for  predicting  the  motion  of,  and  other 
quantities  associated  with, a small  cavity  set  into  motion  by  an 
acoustic  pressure  field  in  a liquid  was  derived.  This  formulation 
consists  of  a set  of  differential,  integral  and  algebraic  non-linear 
equations  that  take  into  account  compressibility  and  viscosity  of 
the  liquid,  heat  conduction  across  the  interface  and  surface  tension, 
The  equations  have  been  programmed  for  solution  on  a digital  com- 
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puter.  In  the  second  paper  , hereinafter  referred  to  as  CD: II, 
this  formulation  was  used  to  study  the  free  pulsations  of  small 
argon-filled  cavities  in  water.  The  results  presented  in  CD: II 
(and  in  an  expanded  report'1'8)  are  fundamental  to  the  interpretation 
of  calculations  reported  here. 

The  limitations  on  the  usefulness  of  such  a mathematical  model 
include  the  assumptions  that  the  speed  of  sound  in  the  liquid  is 
constant  (that  is,  the  variational  pressure  is  a linear  function  of 
the  variational  density) , that  the  cavity  retains  its  spherical 
shape  throughout  its  motion,  that  the  amount  of  gas  in  the  cavity 
remains  constant  and  that  the  gas  in  the  cavity  behaves  as  an  ideal 
gas. 

In  the  calculations  reported  here,  a pre-existing  seed  of  argon 
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in  water  grows  into  a cavity  under  tne  influence  of  a uniform, 
sinusoidal  acoustic  field  characterized  by  a pressure  amplitude,  P^, 
and  a frequency,  f 

The  system  of  notation  adopted  in  CD: I will  be  ured  in  this 

paper.  An  asterisk  (*)  is  used  to  denote  a quantity  in  some  con- 

, * 
sistent  set  or  units.  In  this  notation,  p is  a pressure  in  bars 

* * . * 
and  p = p /p  is  the  non-dimensional  pressure  where  p^n  is  the 

reference  pressure  (here  taken  to  be  1 bar) . Thus  Pa  is  the  non- 

'k  k **  k k 

dimensional  pressure  amplitude,  • The  quantity  R = R /R 

is  the  non-dimensional  radius  of  the  cavity,  R its  radius  in 

* 

centimeters  and  R the  initial  radius  of  the  cavity  in  centimeters w 

n * * £ 

The  frequency,  f,  is  the  non-dimensional  frequency  f t where  f 

k k k k “ 

is  the  frequency  in  Hertz  and  t = Rn/an  and  an  is  the  equilibrium 
speed  of  sound  in  the  liquid. 

Once  set  into  motion,  the  cavity  passes  through  a series  of 

maxima  and  minima  and  its  motion  may  or  may  not  be  periodic.  The 

complexity  of  this  non-linear  motion  results  from  the  tendency  of 

the  cavity  to  pulsate  both  at  the  driving  frequency,  fA,  and  at 

some  resonance  frequency,  f , of  free  pulsation  as  determined  in 

CD: II.  When  the  motion  is  periodic,  there  is  a least  period,  T^,  in 

which  the  motion  repeats  itself.  This  least  period,  T , must  be  an 

Ll 

integral  multiple  of  both  the  acoustic  period,  T^  = 1/f^,  anc^ 

period  of  free  pulsation,  Tr  = 1/f  ? that  is,  for  a periodic  motion, 

Tl  = m Tr  = n T^  where  n and  m are  integers.  Of  particular  interest 

is  the  case  where  m = 1 and  T_  = T = n T. . The  motion  is  then  said 

L r A 

to  contain  a subharmonic  order  n (or,  in  terms  of  frequency, 

fL  = fr  = f A/n  and  there  is  said  to  be  a subharmonic  of  order  1/n) . 

Most  motions  of  cavities  in  an  acoustic  pressure  field  are 
quasi-periodic;  that  is,  the  radius-time  curve  shows  a slowly 
varying  time  interval  that  approximates  a least  period,  T_ , and  the 
maxima  and  minima  in  such  a period  change  slowly  in  amplitude  and 
phase  from  one  such  quasi-period  to  the  next. 

The  non-linear  motion  of  a cavity  is  thus  a combination  of  a 
free  pulsation  and  a driven  pulsation,  both  of  which  contain  a 
fundamental  and  associated  harmonics.  The  driven  pulsation  has 
the  period,  T^,  of  the  acoustic  field,  but  the  period  of  the  free 
pulsation,  T , depends  on  the  amplitude  of  motion,  as  shown  in  CD: II. 
In  general,  in  a least  period,  TL,  the  radius  time  curve  will  be 
quite  complicated.  However,  when  the  amplitude  of 
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motion  is  such  that  the  period,  Tr,  of  free  pulsation 

equals  the  acoustic  period,  T^,  then  the  motion  consists 

of  a single  maximum  and  a single  minimum  in  a least  period, 

T_  = T = T . The  frequency  f at  which  this  coincidence  takes 

place  is  the  non-linear  resonance  frequency  for  a cavity  of  initial 
* 

radius,  Rn,  in  an  acoustic  field  specified  by  the  pair  (PA,  Tft) . 
This  non-linear  resonance  frequency  is  always  less  than  the  linear 
resonance  frequency,  fQ,  for  a cavity  pulsating  with  very  small 
change  in  radius. 

Curves  of  the  resonance  frequency,  fr,  as  a function  of  the 

maximum  radius,  R , are  shown  in  Fig.  2 and  Reference  10  for 

® * 

various  values  of  initial  cavity  radius,  R . When  a cavity  of 

•k  ” 

initial  radius,  R , is  driven  at  an  acoustic  frequency,  f , these 

n A 

curves  tell  us  the  maximum  radius,  R , at  which  the  specified  f 

O Pi 

equals  some  resonance  frequency,  fr,  of  that  cavity.  In  a non- 
linear pulsation  at  a frequency  close  to  some  f , the  radius-time 
curve  has  an  unique  maximum,  Rq,  corresponding  to  the  specified 
pressure  amplitude,  P^.  This  pressure  amplitude  at  which  the 
acoustic  frequency  f.,  equals  a resonance  frequency,  f , will  be 
called  the  resonance  oressure,  Pr. 

We  shall  find  that  the  resonance  pressure,  P , determines  one 
cavitation  threshold.  Another  cavitation  threshold  is  defined 
through  use  of  a function  called  the  dynamical  threshold  radius, 

R. , described  in  CD: II  for  free  pulsations.  When  the  maximum 
*-  * 
radius,  RQ,  of  an  expanding  cavity  of  initial  radius,  Rn,  exceeds 

the  threshold  radius,  the  cavity  is  a transient  cavity.  When  such 

a cavity  collapses,  inertial  forces  in  the  surrounding  liquid 

generate  rapidly  increasing  kinetic  energy  that  is  either  stored 

in  the  compressible  liquid  or  converted  into  internal  energy  of 

the  cavity  contents.  Ultimately,  the  inward  motion  is  halted  by 

the  pressure  in  the  cavity  and  part  of  the  stored  energy  radiated 

as  a shock  wave.  Most  of  the  phenomena  summed  up  as  cavitation 

activity  are  brought  about  by  transient  cavities. 

Determination  of  the  threshold  radius,  Rfc,  requires  partition 

of  the  acceleration  of  the  cavity  interface  into  two  functions:  the 

pressure  function,  PF,  and  the  inertial  function,  IF.  When  a 

cavity  starts  to  contract  from  a maximum,  RQ,  PF  first  decreases, 

passes  through  a minimum  and  then  increases.  IF  is  a function  of 
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the  maximum  radius,  Rq,  at  the  start  of  collapse,  and  the  value  of 
Rq  for  which  IF  intersects  this  minimum  is  the  threshold  radius,  Ft. 

When  a cavity  of  initial  radius  Rn  pulsates  in  an  acoustic 
pressure  field  of  frequency  f^,  its  maximum  radius  increases  when 
the  pressure  amplitude,  P , increases.  The  value  of  P which  causes 

/V  Pi. 

Rq  to  equal  or  exceed  R^  will  be  called  the  threshold  pressure,  Pt* 

It  is  qualitatively  obvious  that  cavitation  activity  must 
increase  with  acoustic  pressure,  but  must  eventually  decrease.  As 
the  pressure  amplitude  increases  the  average  volume  of  a cavity 
becomes  much  larger  than  its  equilibrium  volume  and  the  cavity 
spends  most  of  a period,  Tt , in  such  an  expanded  state.  The  liquid 
then  becomes  much  more  compressible  and  this  increase  in  compressi- 
bility strongly  moderates  the  violence  of  collapse.  The  exciting 
sound  beam  is  both  scattered  and  absorbed  by  the  increased  cross- 
sections  of  the  cavities  and  the  radiated  shock  waves  from  collapsing 
cavities  will  likewise  be  scattered  and  absorbed  by  surrounding 

cavities.  Because  there  may  be  an  enormous  number  of  cavitation 
3 

events  per  cm.  in  a cavitation  zone,  any  increase  in  the  average 

size  of  cavities  may  have  a drastic  effect  on  cavitation  activity. 

11  6 
Sirotyuk  estimates  that  there  may  be  as  many  as  10  cavitation 

events  per  cm. 3 

The  mathematical  model  of  CD: I used  in  carrying  out  the  calcula- 
tions reported  here  predict  the  motion  of  a single  cavity  in  an 
infinite,  homogeneous  liquid.  With  this  restriction  in  mind,  two 
quantities  have  been  chosen  for  calculation 

they  might  give  insight  into  the  thresholds  observed  in  zones 
containing  many  bubbles  with  a wide  distribution  of  initial  radii. 
These  quantities  are  the  maximum  pressure,  p , in  a collapsing 
cavity  and  the  work,  W£,  done  by  a cavity  on  the  surrounding  liquid 
in  expanding  from  its  minimum  radius,  R , to  a maximum  radius,  R , 

ic  jc  1c  HI  — i O 

This  work  is  W_  = W„/W„  where  W =0.88  kilojoules  mol 
ej  £j  n n 

The  maximum  pressure,  p , determines  the  initial  strength  of 
the  radiated  shock  from  a collapsed  cavity  and  W£  measures  the 
transfer  of  energy  to  the  liquid  by  the  compressed  gas  in  the  cavity. 
In  any  expansion,  most  of  the  work,  WE,  is  done  in  the  initial  stage 
when  R is  close  to  R . Thus,  while  p determines  the  strength 
of  the  shock  front,  WE  determines  the  width  and  magnitude  of  the 
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shock  wave  behind  the  front.  Because  the  least  period  TL,  may  be 
many  times  the  acoustic  period,  T.  , the  work,  W„,  is  here  defined 
as  the  average  work  done  by  an  expanding  cavity  in  an  acoustic 
period,  T^,  the  average  being  taken  over  values  in  the  least 
period,  Tt . 

^ £ 

Fig.  1 shows  the  maximum  pressure,  p , predicted  for  a cavity 

. . . . * —4 

of  initial  radius  R = 5 x 10  cm.  as  a function  of  the  acoustic 

n * 

pressure  amplitude,  P^,  for  three  frequencies  of  the  acoustic 

field.  These  frequencies,  600  kHz,  300  kHz,  and  100  kHz,  are 

approximately  equal  to  f , f /2  and  fQ/6  where  f is  the  linear 

resonance  frequency  of  the  cavity. 

At  the  calculated  threshold  pressure,  Pfc,  there  is  a marked 

change  in  the  maximum  pressure  at  300  kHz  and  100  kHz.  At  600  kHz 

the  inertial  function  IF  always  lies  above  the  minimum  in  PF  and 

P.  is  undetermined, 
t 

We  can  draw  two  conclusions  from  the  location  of  the  resonance 
pressure,  Pr,  on  these  curves.  Most  points  on  these  curves  are 
accompanied  by  an  integer.  This  integer  indicates  the  least  period 
in  terms  of  T- . Thus  n = 1 means  Tr  = T* , while  n = 4 means  that 
Tl  = 4 Ta  and  a subharmonic  of  order  1/4  exists.  The  location  of 
Pr  divides  each  curve  into  two  parts.  Below  Pr  only  one  subharmonic 
could  be  found,  while  above  Pr  there  exists  a profusion  of  sub- 
harmonics of  various  orders.  At  300  kHz  and  100  kHz,  the  curves  of 
* 

Pm  abruptly  decrease  in  slope  in  the  vicinity  of  P , but  at  600  kHz 
the  change  in  slope  is  much  less  marked. 

* 

Fig.  3 shows  the  maximum  pressure,  p , as  a function  of  the 

4c  III 

acoustic  pressure  amplitude,  ?A,  for  three  different  cavities. 

Each  cavity  is  driven  at  a frequency  equal  approximately  to  f /2 

* o 

corresponding  to  its  initial  radius,  Rn.  On  all  three  curves  the 

threshold  pressure,  P , marks  an  abrupt  change  in  the  slope  of  the 

curve  and  in  the  vicinity  of  the  resonance  pressure,  p , there  is 

* 

an  even  more  pronounced  decrease  in  the  slope  of  p as  a function 
* m 
of  P^.  Again,  there  are  subharmonics  in  abundance  above  P^  but 

only  one  below  it. 

The  maximum  pressure  curves  would  lead  us  to  identify  Pfc 
with  the  first  cavitation  pressure  threshold  of  Monakhov  et  al 
and  the  resonance  pressure,  P , at  which  f = f with  the  second 

A 37 

cavitation  pressure  threshold.  At  P the  maximum  radius,  R , is 

37  O 
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the  resonance  radius  for  the  driving  frequency,  fA«  These  con- 
clusions would  appear  to  apply  only  to  cavities  driven  at  frequencies 
well  below  the  linear  resonance  frequency,  fQ. 

The  work,  WE,  per  TA  done  by  an  expanding  cavity  is  shown 
in  Fig.  4 as  a function  of  the  acoustic  pressure  amplitude,  P , 
for  a cavity  of  initial  radius,  R = 5 x 10*"  cm.  For  this  cal- 
culated quantity  it  is  even  clearer  that  Pfc  and  P are  the  first 
and  second  pressure  thresholds,  at  least  for  fA  much  less  than  fQ. 

Fig.  5 shows  W^,  as  a function  of  P^  for  three  different 
cavities.  Each  cavity  is  driven  at  a frequency  fft  approximately 
equal  to  fQ/2.  For  each  c&vitv,  p^  is  the  pressure  threshold  at 

which  there  is  a marked  decrease  in  the  slope  of  the  curves.  For 
-5 

the  5 x 10  cm.  cavity,  P.  does  not  appear  to  act  as  a threshold 

* z 

while  the  point  at  P^  = 4 bars  does.  The  significance  of  this 
remark  lies  in  the  fact  that,  at  this  pressure,  a quantity  called 
the  energy  dissipation  modulus,  AW/Wm,  defined  in  CD: II,  is  a 
maximum.  Again,  in  the  curves  for  W£,  subharmonics  appear  in 
general  only  above  the  second  threshold,  Pr. 

The  calculations  reported  here  predict  the  behavior  of  a 
single  bubble  in  an  infinite,  homogeneous  liquid,  and  one  must 
be  cautious  in  seeking  quantitative  correspondences  with  experi- 
mental results,  which  in  general  are  statistical  averages  over 
many  bubbles.  The  maxima  characteristic  of  cavitation  activity  do 
not  appear,  nor  should  we  expect  them  to.  However,  the  results 
clearly  give  us  useful  insights  into  experiments  such  as  those 
of  Monakhov  et  al.  Thus,  despite  these  restrictions,  there  are 
several  general  remarks  that  can  be  made  about  the  results  reported 
in  this  paper: 

1.  The  quantity,  Pfc,  is  a pressure  threshold  at  which  cavitation 
actively  rapidly  increases  for  any  driving  frequency,  fA,  well 
below  fQ,  the  linear  resonance  frequency  of  a cavity.  Tentatively, 

Pt  may  be  identified  as  the  first  cavitation  threshold  of  Monakhov 
et  al. 

2.  On  the  other  hand,  P,  does  not  appear  to  be  a cavitation 

*"  * 

threshold  when  f_  is  approximately  equal  to  fn  or  R is  less 

a ^ u n 

than  a micron.  When  R is  less  than  a micron,  the  first  cavitation 

n 
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threshold  may  occur  when  the  energy  dissipation  modulus,  AW/Wm 
is  a maximum,  as  suggested  in  CD: II. 

3.  A second  pressure  threshold,  P , occurs  when  the  driving 
frequency,  fa,  equals  a non-linear  resonance  frequency,  f , of  the 
cavity.  Tentatively,  P may  be  identified  with  the  second  threshold 
of  Monakhov  et  al. 

4.  At  the  second  threshold,  Pr,  both  the  maximum  pressure  and  WE 
change  abruptly  for  driving  frequencies  less  than  the  linear 
resonance  frequency  of  the  cavity.  Curves  of  both  maximum  pressure 
and  WE  tend  to  flatten  out  for  pressure  amplitudes  greater  than  ?r. 
Changes  in  the  medium  due  to  expanded  cavitation  bubbles,  noted 
above,  may  cause  these  curves  to  decrease  above  P . 

5.  Subharmonics  in  general  are  present  only  in  the  region  above 
the  second  threshold,  P^,  which  may  be  identified  as  the  region  of 
reduced  cavitation  activity  defined  by  Ilonakhov  et  al. 

6.  Cavities  with  initial  radii  greater  than  a micron  are  more 
effective  in  producing  cavitation  activity  when  driven  at  frequencies 
less  than  their  linear  resonance  frequencies. 
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FIGURES 

Fig.  1 Maximum  pressure  in  a 5-micron  cavity 

Fig.  2 Non-linear  resonance  frequency  curve  for  a 5-micron  cavity 

Fig.  3 Maximum  pressure  in  three  cavities 

Fig.  4 Work,  WE,  done  by  an  expanding  5-micron  cavity 

Fig.  S Work,  W^,  done  by  three  expanding  cavities 


13 


Fig.  2 Non-linear  resonance  frequency  curve  for  a 5-micron  cavity 


( bars ) 
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WE,  WORK  DONE  BY  EXPANDING  CAVITY  PER  TA 
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4.  Rn=5xl0"3  cm. 
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Fig.  5 Work,  WR,  done  by  three  expanding  cavities 
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